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coefficient  magnitude,  and  depicting  shock  wave  locations.  On 
the  basis  of  these  results,  it  can  be  verified  that  our 
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I .  INTRODUCTION 

The  flight  of  an  object  over  a  wide  range  of  speeds  has  a 
critical  transition  zone  where  both  subsonic  and  supersonic 
types  of  flow  exist.  This  speed  regime  is  referred  to  as  the 
transonic  range.  The  critical  aerodynamic  behavior  occurs  in 
the  range  0.8  <  M  <  1.2  depending  on  the  object,  where  the 
aerodynamic  coefficients  have  been  found  to  change  by  as  much 
as  100%.  The  behavior  of  the  aerodynamic  force  components  is 
usually  characterized  by  a  rapid  increase  in  the  coefficients 
followed  by  a  sharp  drop,  in  other  words,  a  peak  value  for  the 
coefficients  arises  in  this  regime  [Ref.l  and  2]. 

One  of  the  aerodynamic  components,  aerodynamic  drag, 
represents  a  significant  adverse  force  on  all  flying  objects 
such  as  aircraft,  missiles  and  projectiles.  A  high  drag  force 
reduces  the  craft's  range  capability  or  equivalently  requires 
more  energy  to  achieve  a  certain  range.  Any  effort  to  reduce 
the  drag  coefficient  in  the  design  process  must  concentrate  on 
reducing  the  wake  and  pressure  drag  (inclusive  of  wave  drag) 
contributions  to  the  total  drag. 

A  traditional  approach  to  investigate  aerodynamic 
characteristics  is  based  on  wind  tunnel  data  and  actual  flight 
testing  as  a  source  for  improvement  of  configurations  to  get 
the  best  results.  Unfortunately,  both  wind  tunnel  and  flight 
tests   are   considerably   expensive   and   time   consuming, 


particularly  in  the  transonic  flow  regime.  Because  of  the 
great  lateral  influence  of  profiles  at  transonic  speeds  (shock 
reflection  of  the  wall  bounces  back  to  the  model)  ,  models  must 
be  made  extremely  small  compared  with  wind  tunnel  dimensions 
and  this  introduces  great  experimental  difficulties. 
Consequently,  the  data  obtained  at  transonic  speeds  is 
considerably  less  reliable  than  at  either  subsonic  or 
supersonic  speeds.  In  addition,  from  a  mathematical  point  of 
view,  even  the  two  dimensional  small  perturbation  potential 
equation  for  transonic  flow  retains  one  non- linear  term  which 
is  essential  for  non-divergent  solutions  at  Mach  one,  but  this 
non- linear  partial  differential  equation  has  proven  to  be 
difficult  to  solve.  Such  inherent  difficulties,  coupled  with 
the  presence  of  shocks  in  the  flow  which  cause  boundary  layer 
separation,  have  resulted  in  the  creation  of  many  approximate 
methods  of  solution  which  are  employed  in  the  design  of 
transonic  airfoils  and  the  like  [Ref  .3,  4  and  5]  .  On  the  other 
hand,  the  use  of  numerical  simulation  known  as  Computational 
Fluid  Dynamics  (CFD)  to  predict  aerodynamic  characteristics 
greatly  increases  possibilities  to  improve  design  optimization 
at  relatively  low  cost  and  allows  for  ease  of  design  changes. 
Finally,  using  the  latest  capabilities  of  Euler  as  well  as 
Navier- Stokes  f low- solvers ,  it  has  been  possible  to  compute 
the  flow  over  axi- symmetric  bodies  with  various  contours  in 
the  transonic  regime. 


II  TRANSONIC  FLOW 

Transonic  flows  are  characterized  by  the  simultaneous 
presence  within  the  flow  field  of  both  subsonic  and  supersonic 
regions.  The  properties  of  transonic  flows  can  be  derived  from 
the  general  equations  of  gas  dynamics,  namely,  the  equations 
of  state,  continuity,  momentum  and  energy.  The  following 
derivation  based  on  assumptions  that  the  flow  is  steady, 
irrotational  and  isentropic  with  no  energy  transfer,  no  body 
forces  and  no  shear  stresses  (inviscid  flow) . 

1.   Small  Perturbation  Theory- 
Starting   with   the   equation   of   motion   for   steady, 
isentropic,  inviscid  flow  in  the  index  notation  form  [Ref .4]  : 


Small  -  Perturbation  Theory  gives  the  velocity  field  as 

u^   =  U  +  u       ;   u^  =  V      ;   u^   =  w  (2) 

where  U  denotes  the  free  stream  velocity  in  the  x  direction 
and  u,  V,  w  are  called  perturbation  velocity  components  in  the 
direction  x,  y,  z  respectively. 


Writing  eq.  (1)  out  in  full,  and  substituting  the  velocity 
field  eq.(2),  gives  the  equation  in  terms  of  perturbation 
velocities . 


\  dx        dy        dzl  dx  dy  dz 

,  rr      \    J  du    ^    dv\    ^        I  dv         dw\    ,        /  „,     .  I  dw         du 


(3) 


From  the  energy  equation  for  a  perfect  gas,  the  speed  of  sound 
(a)  can  be  expressed  in  terms  of  the  perturbation  velocities. 


jU+u)^  +  v^  +  w^    ^   a^       =  i^  +  ^' 


2 


(Y-1)     2     (Y-1) 


(4) 


or 


a2  =  ai  -  AUd.  {2uU  ^  u^  +  v^  +  w^)  (5) 


Substituting  eq.(5)  into  eq.(3),  dividing  by  aj^ ,  and 
rearranging  the  terms,  gives  the  full  exact  equation  in  terms 
of  perturbation  velocities  and  free  stream  Mach  number.  This 
equation  contains  linear  terms  on  the  left-hand  side,  but  on 
the  right-hand  side  the  terms  are  nonlinear. 


ox        dy         oz 


m: 


(Y^i)^  +  iYiilii!  ^   (Y-i: 


[/ 


[/' 


in 
ax 


+  M' 


(Y-l)-^    +   jYlllj^l    +     (y-1)     iw^  +  U' 


U 


W 


U' 


dv 
dy 


(6) 


+  m: 


-l)_y   +     (y-^l)    i^^    +     (y-l)     (u^  +  v^) 


(Y-l) 


f/ 


C7^ 


f/^ 


dw 
dz 


+  ML 


!:[      Ij^dy     dxj         U\       ul\dz    'dxj  ^     uA~^    ~dz 


If  the  perturbation  velocities  are  small  (u/U,v/U,w/U  <<  1), 
eq.(6)  can  be  simplified  by  neglecting  the  terms  containing 
squares  of  the  perturbation  velocities  on  the  right  hand  side, 
yielding  : 


ox        dy         dz 


Mi(Y-l)4t|^  .Mi(Y-l)^/|^-^ 
U  ax  U\  dy     dz 

"  U\  dy     dxj         "  U\  dz     dx 


(7) 


For  further  simplification,  in  eq.  (7)  all  the  terms  on  the 
right  hand  side  can  be  neglected,  in  comparison  to  those  on 
the  left  hand  side.  This  gives  the  linear  equation,  which 
contains  only  perturbation  velocities  and  is  valid  only  for 
subsonic  and  supersonic  flow. 


r2\  du        dv    .  dw 


{1-Mt)^    ^  ^  ^  ^  =  0  (8) 

ox        dy         oz 


For  transonic  flow,  where  Mj„  ^  1 ,  the  coefficient  of  du/dx   on 
the  left  hand  side  becomes  very  small,  but  it  is  not  correct 
to  neglect  the  first  term  on  the  right  hand  side  of  eq. (7) . 
Therefore,  the  governing  equation  of  transonic  flow  in  term  of 
perturbation  velocities  is  as  follows: 


,-  ,jr2.   du    ^    dv  ^    dw       .,2  I    ^.  .    u  du  #«» 

(l-M„)-jj_  +     +     =  M„  (7  +  1)  (9) 

ax        ay        dz  U  ax 


For  irrotational  flow,  a  perturbation  velocity  potential  0 
exists, 

u  =  -^  ,  v=  -^  ,  w=  ^                              (10) 

ax  ay                   az 

Substitution  of  eq.  (10)  into  eq.  (9)  ,  gives  the  governing 

equation   for   transonic  flow   in   terms   of   the  velocity 
potential . 


r2 


1  m')-^  -^  -^  -  i!^  -  ^-^'^"^^  d^d^        (11) 

"   3x2     3^2     3^2         t7     ^a^2 


For  bodies  of  revolution,  it  is  convenient  to  use  cylindrical 
coordinates  {x,r,d)  where  x  is  aligned  with  the  body  axis.  The 
velocity  components  corresponding  to  (x,r,^)  are  u^^,  u^,  and 
Uq,  respectively.  The  velocity  potential  in  the  cylindrical 
coordinates , 

".=  ^-=||.   ".=  ||.   u,-l^  (12, 

Transforming  Cartesian  coordinate  eq.(9)  into  cylindrical 
coordinate,  gives  the  governing  equation  of  transonic  flow. 


(1  M^)i!^  *  i^  -  l  ^^  -  J^l!^  =  ^"^^^^^  30 _a^  (13) 

"  3x2    5^2    r-^        ^2  QQ2  u        '^dx^ 


For  axially  symmetric  flow,  where  the  conditions  are  the  same 
in  every  meridian  plane,  there  is  no  variation  with  6,    so  the 
small  perturbation,  non- linear,  axi- symmetric  transonic 
potential  equation  can  be  written  as  follows  : 


d^    ^    Idcp    ^      ^_^2      d^    _    M^iy-1)    d<t>d^  (14) 

3^2    r  "5?        "  a^2        u         ^dx^ 


Or  rewriting  in  shorthand  notation, 


2 


Reference  6  introduces  a  modified  potential  equation  for  axi' 
symmetric  flows 

l-^(r0^)  -  n-M^)4>^  =  <t>J>^  (16) 


where  the  modified  velocity  potentials  are  : 

0^  ^  M^  (7  +  1)  ^   ,         4>r^Ml  (7  +  1)  ^ 


(17) 


2.   An  exact  solution  for  axi- symmetric,  transonic  flow. 

Solutions  to  the  modified  transonic  eq.  (16)  have  been 
given  by  Biblarz  [Ref.6  and  7]  by  using  the  separation  of 
variables  approach  with  a  potential  function  0(x,r)  of  the 
form 

4>{x,r)    =^{x)r]{r)    +  il-M^)  x  (18) 


Substituting  the  above  0  function  in  the  modified  transonic 
eq.(16)  results  in  two  ordinary,  second  order,  non-linear 
differential  equations 

8 


^^  -X?  =0 


dx 


(19) 


dx' 


and 


^  +  1^  -  Xr,2  = 


dr' 


r  dr 


(20) 


where  X  is  the  separation  constant. 

The  solution  to  the  first  differential  eq. (19)  is  obtained 
by  multiplying  both  sides  d^/dx, 


or 


dx 


dld^ 
dx  (^2 


dx 


(XO  =  0 


(21) 


_d_ 
dx 


3  dx 


2^ 


=  0 


(22) 


Thus 


dx     2  ^ 


(23) 


where  o;  is  a  constant  given  by  Ref .  6 
Rearranging  eq.(23) 


dx  = 


d^ 


2    ^ 


Oi 


(24) 


and  integrating  eq. (24)  gives, 


X  -  Xq  = 


d? 


Im^  -M' 


(25) 


where  oi   and  Xq  are  integration  constants. 

The  solution  to  the  second  non- linear  ordinary- 
differential  equation  (20) ,  is  obtained  by  an  outer  expansion 
method.  [Ref.8] 


r]{r)    =  -A-   +  (l-Mi)fi(r)  +  {i-M^)^f^{r)    + 


(26) 


Where  (1-M„^)  represents  a  small  parameter  and  the  first  term 
is  the  purely  sonic  solution. 

By  taking  first  and  second  derivatives  eq. (26)   and 
subtituting  them  into  eq. (20)  yields, 


Vir)  =■ 


\r^ 


28+8\/8" 


(27) 


Al-Hashel  [Ref.9]  reported  on  eq. (26)  and  eq. (27) , 
implemented  the  boundary  conditions  with  the  constants  ol,  C-^, 
a  and  X,  gave  the  final  results  in  the  new  variables. 


?  -  ^ 


3X 

2Ci 


(28) 


10 


_^(yF+2)  _  ^(yF+2)  ^  ^ 


(29) 


X  =  C-, 


1       1 
"6  J  3X\3 


# 


(30) 


>7  = 


V   X 


(a  X 


0.4142 


(31) 


then 


X  = 


d^ 


^2+  \1-M, 


2  I  1.7574 


(32) 


where  Xq=0,  and 


V   = 


4^  |1-mJ|  r(^^^2) 


:i-M^)^ 

50.63 


^2  (yF+2) 


(33) 


Equation  (32)  has  been  numerically  integrated  and  plotted  on 
figure  (1)  as  |  (x)  versus  x,  and  eq.  (33)  will  be  evaluated  and 
plotted  in  figure  (2)  as  rj  (r)  versus  f  for  M^  =  1.05,  1.1  and 
1.2.  A  "patching"  technique  discussed  in  Ref.9  has  been  used 
here. 
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III.  PRESSURE  COEFFICIENT  AND  BOUNDARY  SURFACES 

1.  Pressure  Coefficient 

Liepmann  and  Roshko  [Ref .4]  define  the  pressure  coefficient 
as , 


^P  = 


ip-p^ 


0.5p„C72    yM^ 


^-1 

Poo 


(34) 


From  the  isentropic  relation,  we  have  pressure  ratio  in  terms 
of  Mach  number  and  after  substituting  into  eq. (34)  yields 


c„  = 

2 
7W« 

• 

2+(7-l)M^ 

^p 

2+{y-l)M^ 

Y/(Y-i: 


-  1 


(35) 


Introducing  M^^^  =  Xp" 1 3.^ ,  M^  =  u^/a^  and  using  energy  equation 
(recall  eq.4  and  eq.5),  the  pressure  coefficient  can  be 
expressed  in  terms  of  the  perturbation  velocities. 


^P  = 


1^1 


L.ni 


m: 


1-. 


17+ u)  ^+v^+w^ 
C72 


1y/(y-i) 


1^1 


7-ln.r2 


2u  ^  u^+v^+w^ 


lY/(Y-i: 


-  1 


(36) 
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Using  the  binomial  expansion  on  the  expression  inside  the 
square  brackets  in  eq.  (36)  ,  we  obtain  the  pressure  coefficient 
in  the  form 


^P  =  - 


(37) 


For  axi- symmetric  flow,  in  cylindrical  coordinates  where  u=Uj^ 
and  (v^+w^)  =  u^.^  ,  substituting  into  eq.  (37)  yields 


^P  =  - 


U„ 


(38) 


The  linearized  pressure  coefficient  approximation  for  axi 
symmetric  flow  turns  out  to  be 


C    = 


(39) 


Recall  the  modified  axial  velocity  potential,  eq. (17) , 


0^  =  ^^('y  +  l)^ 


(40) 


thus 


4>: 
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(41) 


Substitute  eq.  (41)  into  eq.  (39) ,  yields 


P       "72 


^00(7^1 


(42) 


The  derivative  of  the  potential  function  eq. (18)  with  respect 
to  X  becomes 


0^  =  ^  ^  -  (i-^c^: 


(43) 


Ref.6  introduces  relation  of  the  constants  C^^,  01,    a  and  X  in 
expression, 


q;  =  tCjl-Mfji-^^-^^ 


(44) 


and 


-,1.2426  _ 

^      =  1.08x10  "2 


0.7574 


(45) 


Rewriting  eq.  (23)   and  substitute  the  constant  o;  eq.  (44) 
becomes 


^  - 
dx 


^e  ^   Cil  (1-M, 


l)\^-'^'^ 


(46) 
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Recall  the  new  variable  eq. (28 


^  =  ^ 


3X 


^  1 


2Cn 


(47) 


Substituting  eq.  (47)  into  eq.  (46)  and  factoring  C^^  out,  yields 


dx 


C. 


r  -  I (i-M^) I 


2n  I  1 . 7574 


(48) 


Recall  and  arrange  new  variable  eq. (3i; 


(a  X)0-4i42  ^ 

^  = X ' 


(49) 


Rewriting      eq.(43)       the      modified      potential      function      and 
inserting   eq. (48)    and  eq. (49) ,    becomes 


'aX 


0  .4142 


4>^=   1^^ rj   C,^    [^^.\a-M:)\ 


r2v    I  1.7574 


1-M 


2,        (50) 


or 


^0.4142^.^3 
0x  =    ^-^r-FTT^i^    V 


X 


0.5858 


rH  (iX)l 


2n    I  1.7574 


l-Ml 


(51) 
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Substituting   the  expression  of   constants   eq.  (45)    into  eq.  (51! 


(p^  =  0.2208  rj 


|2    +    I  il-M, 


2n    I  1  .7574 


^  ^  ii-Mh 


(52) 


Finally,  rewriting  the  pressure  coefficient  C_,  as 


C. 


-2 


P  — 2 


3.220arj 


I'M, 


2x    I  1.7574 


^  +  (i-Mi 


(53) 


2 .   Boundary  Surfaces 

For  an  inviscid  flow,  the  condition  to  be  applied  at  the 
surface  of  a  solid  boundary  is  that  the  direction  of  the  flow 
velocity  vector  must  be  tangent  to  the  solid  surface.  In  other 
words,  the  velocity  vector  is  everywhere  at  right  angles  to 
the  normal  to  the  solid  boundary  [Ref  .10]  .  In  addition  the 
boundary  condition  requires  that  the  gradient  of  potential  (p 
vanish  far  ahead  of  the  body.  In  terms  of  perturbation 
velocities  this  boundary  condition  becomes 


dr 
dx 


suzface 


(54) 
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Recall  the  modified  velocity  perturbation  potential  eq. (17; 


4>r  =  mI    (7^1) -j/  (55) 


thus 


^«    mI    (7  +  1) 


substituting  eq. (56)  into  eq. (54) , 


dr 


(56) 


^'   Wf  (7-1) 


(57) 


Taking  the  differential  of  eq. (18)  with  respect  to  r  gives 


0r  =  ^  -g  (58) 


Expressing  eq. (58)  in  terms  of  new  variables,  recall  eq. (28) , 
eq. (29)  and  eq. (31) 


4>r  =   0.0848  f  -^  (59) 

dr 
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Recalling  eq.  (33)  and  taking  the  derivative  with  respect  to   r  , 
then  substituting  eq.(59)  into  eq. (57) ,  yields 


dx 


0.0848 

Mi(Y+l) 


I- 


+  2.8284     \1-MI\     f^-^284+Q      ]^5^2  (1-Mf)2     f6.657 


(60) 


Replacing  the  left  hand  side  of  eq. (60)  with  the  new  variables 
eq. (29)  and  eq(30),  then  rearranging  it  becomes 


di 
dx 


0.0326 
Mi(Y+l) 


^+2.8284  \1-mI\   f^-^284^.Q_i5i2(i-Mi)2  ^6.657 


(61) 


For  further  arrangement,  by  separating  variables  eq.  (61) , 
yields 


df 


— -2.8284|l-Mi|  fi-8284_o.i5i2(l-Mi)2  ^e.esv 


f^ 


(62) 


-0.0326 
Mi(Y+l) 


I   dx 


Then,  integrating  both  sides  of  eq. (62)  becomes 
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/ 


dr 


-2.8284|l-Mi|f  ^-^284  -0.1512  (1-MJ)2  f^-^^^ 


0.0326  C  r 


-0.0326 


dx 


(63) 


where 


r„  = 


1.2074 

|1_^|0.2071 


(64) 


Al-hashel  [Ref.lO]  has  developed  and  computed  eq.(63) 
using  numerical  integration,  to  determine  the  boundary- 
surfaces  in  dimensional  and  non-dimensional  (normalized)  form 
for  M„  =  1.05,  1.10  and  1.20  as  depicted  in  Figure  (3)  and 
Figure  (4) .  Based  on  these  calculations,  the  geometric  grid  of 
the  afterbody  for  Mach  number  1.10  and  1.20  are  developed  for 
further  study  with  Computational  Fluid  Dynamic  (CFD) . 

This  thesis  research  also  examines  conical  afterbodies  as 
a  solid  afterbody  boundary  surface  with  base  diameter  ratio 
(dj^/dj^)  of  0.50  and  0.75  and  conical  turning  angle  (/5)  of  26.6 
and  14.0  degree  respectively  [Ref.ll  and  12].  Then,  for 
further  investigation,  the  complete  bodies  as  a  solid  boundary 
surface  are  generated,  with  a  kind  of  forebody  (ellipsoid) 
joined  with  the  contour  surface  afterbodies  as  well  as  the 
conical  afterbodies. 
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IV.   DRAG 

Drag  is  one  of  the  aerodynamic  force  components  parallel 
to  free  stream  velocity  {U^)  .  It  represents  a  significant 
adverse  force  on  all  flying  objects.  Basically,  the  drag  force 
is  divided  into  two  categories,  the  drag  caused  by  forces 
acting  normal  to  the  boundary  surface  which  is  called  pressure 
drag  (inclusive  forebody,  base  and  wave  drag)  and  that  arising 
from  the  tangential  forces  acting  on  the  surface,  by  virtue  of 
viscosity  which  is  called  viscous  drag  or  skin  friction. 
Figure  (5)  shows  the  components  of  the  drag  as  function  of 
Mach  number  and  the  methods  used  to  compute  body  drag  in  four 
Mach  number  regimes  [Ref .5] . 


Fie 
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b#r     of    .IrH    componrnt*. 
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Fif        >       B^Hv    m'»tho.l«     for    pr 


r»    <l»«f 


Figure  5.   Components  of  the  Drag  [Ref. 5] 
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1.   Skin  Friction 

Skin  friction  is  the  result  of  shearing  force  within  the 
boundary  layer  of  a  viscous  flow,  acting  tangentially  to  a 
surface  in  motion  relative  to  the  fluid.  The  amount  of  viscous 
resistance  depends  upon  whether  the  flow  is  laminar  or 
turbulent  [Ref.5].  Krasnov  [Ref.l3]  introduces  the  laminar  and 
turbulent  skin  friction  for  flat  plates,  based  on  the  boundary 
layer  theory,  for  compressible  flow.  These  formulas  also  valid 
for  bodies  of  revolution  with  infinite  length. 


Laminar  flow,  Re  <  10^ 


CfL  =      ^'^^  (1  ^  0.03M^)''  (65) 


(Re) 


2 


Turbulent  flow,  Re  >  10^  : 


Crr  -         "-^''3   [log.oi^e  (1^0 .  2Mi)  — f"^^       (gg) 
(1+0. lAC) 


where  Re  is  the  Reynold  number 
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2 .   Pressure  Drag 

Liepmann  and  Roshko  [Ref.4]  introduce  the  pressure  drag 
formula  (inclusive  wave  drag  and  base  drag)  for  axi- symmetric 
bodies  as, 


D  =  J  p  dS  -  PbS{L) 


=  /  ip-pJdS   +  ip^-Ps)S{L) 


or,  in  dimensionless  form, 


(67) 


D 


^   g«5(L) 


S{L) 


/ 


^  dx  P^ 


(68) 


where  S  (x)  =  7rr^  is  the  cross-sectional  area  of  the  body  at  x, 


dSix)    =  2nrdr  =  2Tzr—dx 

dx 


(69) 


or 


dS   _  „^^  dr 

— -  -  2zr—- 

dx  dx 


(70) 
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The  first  term  in  square  bracket  eq.  (68)  can  be  solved  by- 
inserting  eq.  (53) ,  eq.  (60)  and  eq.  (70)  yielding, 


L  L 

P~dbc 


-2 


M„(7>1) 


0.2208^[p+|  (l-M^)  |i-7574j   3^(i_^2, 


(71) 


(27rr)' 


2jL^mf  t^  ^2  .  8284  I  1-Mf  I  r^ -8284^0  .  1512  (1-mJ)  2r6-657 


M^(7*l)  \: 


f3 


dx 


Recalling   new  variables   eq. (28) ,   eq. (29) ,   eq. (30)   and 
inserting  |3^  =  l-M^^^  into  eq.  (71)  yields, 


-pI^  = 


-2 


(l-iS^)  (7+1) 


0.0326 


0.2208i?[|+|i33-5i48|jT^^2 


5.1302Ci 


^0.7071  ^0.2071 


(72) 


^+2.8284|i32|ri-8284^Q^15-L2i3*r6-657 


\dx 


where  C-^,    a  and  X  are  constants 
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The  second  term  in  square  bracket  eq. (68)  is  the  contribution 
to  the  drag  by  the  base  pressure  pg.  Ref.4  claims,  that  the 
values  of  the  base  pressure  coefficient  C  3  must  be  obtained 
experimentally.  However  Krasnov  [Ref.l3]  introduces  the 
boattail  drag  coefficient  for  a  conical  afterbody  as, 


CpQ  =   0.002 


0.8 


Ml 


,1.7 


1-S, 


(73) 


where 


6        =  turning  angle  of  the  conical  body 

Sg   =  ratio  of  the  base  area  to  the  mid- section  area 


In  this  work,  the  skin  friction  coefficient  as  well  as  the 
pressure  drag  coefficient  of  the  geometries  of  interest,  at 
specified  Mach  number,  can  be  obtained  directly  from  the 
results  of  CFD  calculations.  The  Euler  solution  will  give  the 
pressure  drag  and  the  Navier-Stokes  solution  will  give  both. 
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V.  COMPUTATIONAL  FLUID  DYNAMICS 

The  rapid  advancement  in  the  speed  of  computers  and  their 
enormous  memory  size  has  led  to  the  emergence  of  the  field  of 
computational  fluid  dynamics  (CFD) .  This  branch  of  fluid 
dynamics  complements  its  experimental  and  theoretical 
branches,  by  providing  an  alternate,  cost-effective  means  of 
simulating  real  flow.  It  also  offers  the  only  means  of 
examining  theoretical  advances  for  conditions  unavailable 
experimentally.  As  a  model  based  method,  CFD  can  provide  the 
convenience  of  being  able  to  switch  off  specific  terms  in  the 
governing  equations  [Ref.l4],  so  as  to  assist  the  researcher 
in  understanding  the  contributions  of  various  physical 
factors. 

In  this  work,  CFD  was  used  to  compute  the  axi- symmetric 
flow  over  the  afterbody  geometry  of  models  only  (the  boundary 
surfaces  obtained  by  the  small  perturbation  method  [Ref .9]  and 
conical  afterbody)  and  over  complete  body  models  which  are 
composed  of  forebody  (ellipsoid)  and  afterbody. 

1.   Grid  Generation 

The  computer  programs  GRAPE  [Ref. 15]  and  GRIDGEN2D 
[Ref .16]  are  tools  used  to  generate  two-dimensional  structured 
grids  about  airfoils  and  other  shapes  by  the  use  of  algebraic 
or  Poisson  differential  equation  solvers.  GRAPE  was  used  for 

25 


the  geometry  with  smooth  contour  surfaces,  while  GRIDGEN2D  was 
used  for  the  geometry  with  a  conical  afterbody  with  non- smooth 
points.  Outer  and  inner  boundaries  were  specified  as  the  C- 
type  grid  for  afterbody  models  only,  while  the  0-type  grid  was 
for  the  complete  body  models,  where  both  type  of  grids  treat 
the  surface  of  the  body  as  the  inner  boundary.  The  important 
characteristics  in  a  grid  generation  technique  are  the  ability 
to  specify  the  spacing  between  mesh  points  at  the  boundary,  in 
the  direction  normal  to  the  boundary,  and  the  control  of  the 
angles  with  which  mesh  lines  intersect  the  boundaries  which  is 
known  as  orthogonality  [Ref.15]. 

Figure  (6)  is  a  typical  output  of  program  GRAPE  for  an 

afterbody  with  small  perturbation  solution  contour  for  Mach 

1.10  (SPS_1.1),  C-type  grid  with  the  grid  size  of  115x60,  and 

figure  (7)  is  the  complete  body  with  the  grid  size  of  152x60 

(0-type  grid) . 

To  develop  a  three-dimensional  grid  from  two-dimensional 
grid  (output  program  GRAPE),  the  FORTRAN  code  called  HALF.F 
(Appendix)  was  used  to  write  out  the  half  of  the  2D  grid  and 
convert  them  into  3D  plane -grid.  Then,  the  FORTRAN  code 
ROTATEGR.F  (Appendix)  was  used  to  rotate  the  3D  plane -grid  for 
11  planes  plus  2  more  to  generate  a  3D  volume  grid  of  half  of 
an  axi- symmetric  body  surface  as  shown  in  figure  (8)  and 
figure  (9)  with  grid  size  of  58x13x60  and  77x13x60 
respectively.  These  grids  are  ready  for  further  processing 
with  a  flow- solver. 
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Program  GRIDGEN2D  gives  the  pattern  of  how  to  build  2D- 
grid.  It  consists  of  four  subfaces  which  are  treated  as 
boundaries.  Grid  spacing  is  determined  by  setting  up  the 
distribution  points  of  each  pair  subfaces.  It  uses  equal 
spacing  along  the  body  surface,  while  the  direction  normal  to 
the  body  surface  uses  geometric  spacing  with  specified  width 
in  the  beginning.  Figure  (10)  and  Figure  (11)  are  typical 
outputs  of  GRIDGEN2D  for  an  afterbody  and  for  a  complete  body 
with  the  grid  size  of  58x60  and  77x60  respectively.  FORTRAN 
code  D2D3 . F  (Appendix)  is  used  to  convert  2D  grid  output  from 
GRIDGEN2D  into  3D  plane -grid.  Then  we  apply  the  FORTRAN  code 
ROTATEGR.F  to  generate  a  3D  volume  grid  of  half  of  an  axi- 
symmetric  body  as  depicted  in  Figure  (12)  and  Figure  (13)  with 
the  grid  size  of  58x13x60  and  77x13x60. 

In  this  research,  we  also  attempted  to  develop  a  fine  grid 
for  complete  bodies  (0-type  grid),  where  the  radius  of  outer 
boundary  is  set  up  to  be  five  times  of  the  body's  length.  The 
steps  to  obtain  3D  volume  grid  are  similar  to  those  described 
above  with  the  final  grid  size  of  77x9x120  and  the  final 
geometry  of  a  quarter  of  an  axi- symmetric  body  instead  of  a 
half  body.  Figures  (14)  and  (15)  are  typical  fine  grids  of  the 
complete  body  with  a  conical  afterbody  and  a  complete  body 
with  contour  surface  afterbody. 
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2.   Flow- solver 

The  OVERFLOW  program  [Ref.17]  was  developed  by  NASA  Ames 
Research  Center.  It  uses  either  3-D  Euler  or  Navier- Stokes 
flow- solvers  for  inviscid/viscous  flow,  by  setting  True  or 
False  the  parameter  VISINP  (viscosity  input)  in  the  input  file 
(overflow. in) .  Before  flow-solver  code  (OVERFLOW)  is  applied 
on  the  grid  file,  a  formatted  3D  grid  file  named  "grid. for" 
must  be  converted  into  an  unformatted  grid  file  named 
"grid. in",  by  using  FORTRAN  code  called  READX.F  (Appendix). 
Then,  the  NAMELIST  input  file  parameter  specification  must  be 
written  for  running  OVERFLOW,  it  is  called  "overflow. in" .  The 
input  parameter  consists  of  the  number  of  iterations, 
timesteps,  calculation  methods,  smoothing,  type  of  flow  and 
boundary  conditions  for  each  grid.  The  value  of  angle  of 
attack  (ALPHA)  depends  on  the  orientation  of  the  grid  in  the 
coordinate  system.  In  this  case,  ALPHA  is  180°  (flow  comes 
from  the  x-positive  to  the  x- negative  direction) .  The  boundary 
conditions  depend  on  the  geometry  corresponding  to  the  final 
3D  volume  grid.  Both  input  parameters,  the  angle  of  attack  and 
the  boundary  conditions  for  each  geometric  shape,  are 
tabulated  in  Table  1. 

The  file  overflow,  in  (Appendix)  is  a  typical  input 
parameter  specification  of  the  axi- symmetric  body  with  a 
conical  afterbody  using  the  grid  size  of  77x13x60.  This  file 
input  uses  NSTEP=100,  Mach  number=0.95,  ALPHA=180°.  The 
calculation  method  depends  on  a  central  difference  Euler  term 
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in  J,  K,  L  and  ARC3D  diagonal  factorization  and  dissipation 
scheme  (IRHS=0,  ILHS=2,  IDISS=2).  For  the  first  attempt  a  time 
step  DT=0.1,  ITIME=1  and  CFLMIN=5.0  was  utilized.  The  boundary 
condition  consists  of  6  boundaries  (NBC=6),  with  the  type  of 
BC  IBTYP=30  (inflow)  in  the  J  positive  (1) ,  IBTYP=15  (Axis  K 
round)  in  J  negative  (-1)  direction,  IBTYP=12  (symmetry  in  Y) 
in  the  K  positive  (2)  and  K  negative  (-2)  direction,  IBTYP=1 
(Inviscid  adiabatic  wall)  in  the  L  positive  (3)  direction  and 
IBTYP=32  (Supersonic  /subsonic  inflow/outflow)  in  the  L 
negative  (-3)  direction. 

The  OVERFLOW  program  gives  output  files  such  as  ovr.out, 
q.save,  resid.out  and  fomo.out.  To  verify  that  a  calculation 
is  appropriate  or  converged,  it  can  be  traced  by  looking  at 
the  plot  of  residual  history  (resid.out).  A  convergence 
criterion  was  defined  as  the  reduction  in  residuals  by  two 
orders  of  magnitude.  If  the  first  run  (NSTEP=100)  does  not 
fulfill  the  convergence  criterion,  one  may  perform  a  restart 
by  further  running  OVERFLOW  and  changing  the  parameter 
RESTART=.T.  in  the  input  file  (overflow. in)  and  by  copying 
q.save  into  q, restart.  Repeating  the  above  steps  until  all 
convergence  criteria  are  met.  If  in  the  output  file  q.bomb 
appears,  one  tries  another  run  by  changing  parameter  time  step 
DT  and  CFLMIN  until  the  criteria  are  met.  Then,  one  converts 
the  output  file  q.save  (unformatted)  into  formatted  file  named 
q.form  by  using  FORTRAN  code  READl.F  (Appendix) . 

Finally,  the  last  step  is  to  run  PL0T3D  program  with  the 
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source  formatted  files  x=grid.for  and  q=q.fonn.  By  using 
functions  on  the  PL0T3D  program,  many  plots  such  as  pressure, 
velocity,  Mach  number,  vorticity,  etc.,  are  obtained. 

3.   Results 

The  results  of  CFD  programs  are  grouped  into  the 
corresponding  geometric  shape,  namely,  afterbody  only  and 
complete  body.  The  Euler  flow- solver  (inviscid  flow)  was 
applied  to  all  axi- symmetric  bodies,  except  for  the  afterbody 
models  only,  where  both  Euler  and  Navier-Stokes  (viscous  flow) 
were  applied.  Most  of  the  calculations  converged  in  500 
iterations,  meaning  that  the  residual  history  achieved  a  two 
order  of  magnitude  drop. 

Results  can  be  analyzed  by  plotting  the  q.form  file, 
output  file  from  the  OVERFLOW  program,  using  the  proper 
function  in  PL0T3D  program  (Mach  number) .  By  interpreting  the 
Mach  number  contour  surrounding  the  body  surfaces,  one  can  be 
determine  the  characteristics  of  the  flow  field.  In  addition, 
the  drag  coefficient  C^,  can  be  obtained  in  the  file  fomo.out, 
output  from  OVERFLOW.  Hence,  in  the  sequences  of  Mach  number, 
one  can  describe  the  significant  flow  characteristic  of  each 
geometric  shape. 

a.   Afterbodies 

An  Euler  as  well  as  a  Navier-Stokes  flow- solver  was 
applied  to  these  afterbody  models.  The  approaching  free  stream 
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Mach  number  (ranging  from  1.05  to  1.50)  starts  from  the  mid- 
section of  the  afterbodies,  then  the  flow  follows  along  the 
afterbody  surface  until  it  reaches  a  maximum  local  Mach  number 
as  tabulated  in  Table  2 .  The  maximum  local  Mach  number  for 
viscous  flow  is  always  lower  compared  to  in  the  inviscid  flow. 
This  may  be  caused  by  the  viscous  flow  itself  since  we  are 
taking  into  account  the  shear  force  in  the  boundary  layer  near 
the  surface. 

For  the  afterbodies  from  the  small  perturbation  solution 
contour  (SPSl.O,  SPSl.l  and  SPS1.2),  in  inviscid  flow,  the 
shocks  are  formed  at  the  contour  surface.  The  location  of  the 
shock  depends  on  the  specific  afterbody  contour  and  the 
approaching  Mach  number;  at  the  higher  Mach  numbers  the  shock 
appears  a  bit  further  downstream  as  shown  in  Figures  (16) 
through  (21)  .  As  shown  in  Figure  (22)  and  Figure  (23)  ,  this  is 
apparently  a  result  of  viscosity;  it  shows  the  boundary  layer 
by  the  increment  of  Mach  number  away  form  the  surface.  Weak 
shocks  were  formed  further  downstream  compared  to  inviscid 
flow.  Flow  separation  occurs  in  the  starting  contour  region 
and  is  followed  by  circulating  flow  in  the  base  region. 

For  the  conical  afterbody,  the  approaching  free  stream 
Mach  number  increases  following  the  mid- section  surface,  then 
a  flow  expansion  occurs  at  the  turning  angle  region,  until  the 
maximum  local  Mach  number  is  reached.  A  weak  shock  is  formed 
at  the  end  of  the  boattail  region  as  shown  in  Figure  (24)  , 
while  for  the  viscous  flow,  the  weak  shock  develops  away  from 
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boattail  surface  due  to  the  boundary  layer  and  the  separation 
of  the  flow.  The  circulating  flow  in  the  base  region  is  more 
significant  than  in  the  inviscid  flow  as  depicted  in  Figure 
(25) .  In  addition,  from  Table  2,  for  viscous  flow  the  maximum 
local  Mach  number  for  the  conical  afterbody  is  higher  than  for 
the  contoured  afterbody. 

For  each  afterbody  the  pressure  drag  coefficient  (C^) 
versus  free  stream  Mach  number  (M^^)  for  inviscid  and  viscous 
flow  are  plotted  in  Figure  (26) .  The  negative  sign  of  C^  is 
due  to  the  fact  that  the  calculation  of  pressure  starts  from 
mid- section  through  the  base  of  afterbody  and  ignores  the 
forebody  pressure.  These  results  show  that  the  pressure  drag 
coefficient  is  higher  for  viscous  flows  than  for  inviscid 
flows  for  each  given  afterbody.  This  may  be  caused  by  the 
viscosity  effect  and  the  pressure  distribution  difference  in 
the  flow  field.  In  addition,  it  can  be  seen  from  the  chart, 
that  the  C^  for  the  afterbody  with  small  perturbation  solution 
contour  Mach  I.IO  (SPSl.l)  has  the  lowest  C^  values  over  the 
entire  Mach  number  range.  Therefore,  the  SPSl.l  contour  shows 
to  be  relatively  the  best  among  these  afterbodies. 

b.   Complete  Bodies 

Complete  bodies  consist  of  an  arbitrary  forebody 

(ellipsoid)  joined  to  various  afterbodies  such  as  the  small 

perturbation  solution  contours  (SPS_1.1  and  SPS_1.2)  and  a 

conical  afterbody.  The  approaching  free  stream  Mach  number  {M^ 
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ranging  from  0.50  to  1.50)  starts  from  outer  boundary  with 
distance  of  1  times  for  the  coarse  grid  and  5  times  body's 
length  for  the  fine  grid  respectively.  Significant  differences 
between  coarse  and  fine  grid  calculations  can  be  seen  from  the 
characteristic  bow  shock.  In  the  coarse  grid,  the  bow  shock 
hits  the  outer  boundary,  this  causes  the  approaching  free 
stream  Mach  number  not  to  be  the  same  as  in  the  input 
parameter  as  shown  in  Figure  (27) .  While  in  the  fine  grid,  the 
bow  shock  dies  out  before  reaches  the  outer  boundary  as 
depicted  in  Figure  (28) .  So,  further  discussion  is  focused  on 
the  fine  grid  exclusively. 

The  flow  stagnates  on  the  nose  and  then  follows  the  body 
surface  until  it  reaches  a  maximum  local  Mach  number  as 
tabulated  in  Table  3 .  The  critical  Mach  number  for  these 
complete  bodies  is  approximately  at  M„  =  0.70,  where  the 
maximum  local  Mach  number  reaches  unity  at  the  shoulder 
region. 

Figure  (29)  shows  a  typical  high  subsonic  free  stream  Mach 
number  (M„  =  0.85)  flow  over  the  complete  body  SPS_1.1.  The 
flow  stagnates  on  the  nose  tip,  then  flow  is  accelerated 
following  the  forebody  surface  reaching  a  sonic  line  at  the 
mid-way  of  the  forebody  and  forming  a  supersonic  region  at  the 
shoulder.  Then,  the  flow  deccelerates  at  the  mid-section  and 
accelerates  again  until  reaches  a  maximum  local  Mach  number  at 
the  starting  region  of  the  afterbody  and  a  weak  shock  occurs 
in  this  region.  A  supersonic  approaching  free  stream  Mach 
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number,  is  depicted  in  Figures  (30)  and  (31)  for  SPS_1.1  and 
SPS_1.2  respectively.  A  bow  shock  is  obviously  seen  in  front 
of  the  nose.  The  characteristic  of  the  bow  shock  is  more 
inclined  down  stream  for  the  higher  M^,  as  shown  in  Figure 
(32).  A  subsonic  region  is  formed  between  bow  shock  and  the 
nose,  then  the  flow  accelerates  along  the  forebody  surface  up 
to  a  supersonic  region  in  the  mid- section.  The  expansion  flow 
occurs  in  the  starting  contour  region  until  it  reaches  a 
maximum  local  Mach  number.  Then,  a  shock  is  formed  in  the 
contour  region.  Similar  as  in  the  afterbody  only,  the  shock 
location  depends  upon  the  contour  surface  and  M^.  The  shock 
location  for  a  given  contour  is  more  downstream  for  higher  M^^ 
and  at  the  same  M^,  the  shock  location  for  SPS_l.l  is  more 
downstream  than  SPS_1.2.  A  typical  residual  history  of  CFD 
calculation  for  complete  body  with  small  perturbation  solution 
contour  is  shown  in  Figure  (33).  Convergence  is  obtained  at 
about  500  iterations. 

For  the  conical  afterbody,  the  flow  characteristic  is  the 
same  as  the  other  complete  bodies  up  to  the  mid- section 
region.  The  expansion  flow  occurs  at  the  turning  angle,  then 
the  flow  accelerated  along  the  conical  surface  and  weak  shock 
is  formed  at  the  edge  of  base.  Figures  (34)  and  (35)  show  the 
Mach  number  contour  and  corresponding  residual  calculation  for 
conical  afterbody  at  M„  =  1.10. 

The  pressure  drag  coefficient  (C^^-press)  versus  free 
stream  Mach  number  (M^^)  for  fine  grid  complete  bodies  are 
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plotted  in  Figure  (36)  .  The  drag  rises  sharply  in  the  high 
subsonic  Mach  number  (M^,  =  0.95)  and  reaches  a  maximum  (peak) 
at  M^  =  1.10.  Then  the  drag  decreases  with  a  shallow  curve  as 
the  M^  increases.  The  decreasing  shallow  curve  may  be  caused 
by  the  bluntness  of  the  nose  and  it  agrees  with  Shapiro 
[Ref .18]  because  the  fineness  ratio  and  bluntness  of  the  nose 
of  bodies  of  revolution  are  the  important  factors  that 
contribute  the  drag  curve  at  transonic  and  supersonic  range. 
It  can  be  seen  from  the  graph,  the  drag  curve  of  complete  body 
with  conical  afterbody  is  higher  than  with  small  perturbation 
solution  contour  at  the  entire  M^.  Furthermore,  the  peak  of 
the  drag  curve  is  approximately  15%  higher.  The  drag  curve  for 
SPS_1.1  and  SPS_1.2  are  likely  to  have  the  same  trend  up  to  M„ 
=  0.95;  beyond  this  Mach  number,  the  drag  curve  for  SPS_1.2  is 
slightly  greater  than  for  SPS_1.1.  Therefore,  the  complete 
body  with  small  perturbation  solution  contour  afterbody  Mach 
1.1  (SPS  1.1)  relatively  gives  the  lowest  drag. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  use  of  numerical  simulation  (CFD)  appears  to  be  the 
most  cost  effective  method  to  predict  the  aerodynamic 
performance,  especially  in  the  transonic  range. 

In  this  research,  the  grid- generating  program  GRAPE  is 
suitable  only  for  the  geometry  with  smooth  contour  surface 
(SPS_1.1  and  SPS_1.2) ,  while  the  program  GRIDGEN2D  is  used  for 
the  geometry  of  a  conical  afterbody  with  non- smooth  points. 

We  have  shown  that  for  a  complete  body  model  the  use  of  a 
fine  grid  (77x9x12  0)  is  more  reliable  than  a  coarse  grid 
(77x13x60)  .  This  is  shown  by  the  characteristic  of  the  bow 
shock  at  M(jo  >  1  as  seen  in  Figures  27  and  28.  The  pressure 
drag  coefficient  (C^^)  versus  free  stream  Mach  number  {M^) 
graphs  show  that  the  small  perturbation  solution  contour  for 
Mach  1.10  (SPS_1.1)  gives  relatively  the  lowest  C^  on  both 
models  (afterbody  and  complete  body) ,  a  decrease  of 
approximately  15%  of  peak  in  the  transonic  range  compared  to 
the  conical  afterbody.  Therefore,  the  best  design  of  an  axi- 
symmetric  body  such  as  missiles,  projectiles  and  aircraft,  can 
be  based  on  the  small  perturbation  solution  contour. 

Finally,  this  work  may  be  continued  with  the  investigation 
a  complete  body  with  pointed  nose  geometry  and  the  application 
of  a  Euler  as  well  as  Navier-Stokes  flow-solvers  within 
various  small  angle  of  attack. 

36 


LIST  OF  REFERENCES 


1.  Sahu,  J.,  "Numerical  Computations  of  Transonic  Critical 
Aerodynamic  Behavior,  "  AIAA  Journal,  Vol  2  8  No.  5, May  19  90. 

2.  Sahu,  J.,  "Three -Dimensional  Base-Flow  Calculation  for  a 
Projectile  at  Transonic  Velocity,"  AIAA  Journal,  Vol.  27 
No.  2,  February  19  89. 

3.  Guderley,  K.  G.  ,  "The  Theory  of  Transonic  Flow,"  Pergamon 
Press,  London  1962. 

4.  Liepmann,  H.  W. ,  Roshko,  A.,  "Elements  of  Gas  Dynamics," 
John  Wiley  and  Sons,  Inc.,  1957. 

5.  Nixon,  D. ,  "Transonic  Aerodynamics,"  AIAA,  Progress  in 
Astronautics  and  Aeronautics,  Vol.  81,  1982. 

6.  Biblarz,  0.,  "Phase  Plane  Analysis  of  Transonic  Flows," 
AIAA,  Paper  No.  76-332,  July  1976. 

7.  Biblarz,  0.,  "An  Exact  Solution  to  the  Transonic 
Equation,"  Israel  Journal  of  Technology,  Vol.  13,  1975. 

8.  Van  Dyke,  M.D. ,  "Perturbation  Methods  in  Fluid  Mechanics," 

Academic  Press,  NY  1964. 

9.  Al-hashel,  W.I.,  "Two-Dimensional  Boundary  Surfaces  for 
Axi- Symmetric  External  Transonic  Flows,"  MS  Thesis,  Naval 
Postgraduate  School,  1993. 

10.  Zucrow,  M.  J.,  Hoffman,  J.  D."Gas  Dynamics,"  Vol.  I  &  II, 
John  Wiley  and  Sons,  Inc.,  1976. 

11.  Chow,  W.L.,  "Base  Pressure  of  a  Projectile  Within  the 
Transonic  Flight  Regime,"  AIAA  Journal,  Vol.23,  No. 3, 
March  1985. 

12.  Payne,  P.R.,  "Drag  of  Conical  and  Circular  Arc  Afterbodies 
Without  Jet  Flow, "  Aviation  and  Surface  Effect  Department, 
DTNSRDC/ASED- 80/10,  May  1980. 

13.  Krasnov,N. F. , "Aerodynamics  of  Bodies  of  Revolution," 
American  Elsevier  Publishing  Co., Inc.,  New  York  1970. 

14.  Fletcher,  C.A.J. ,  "Computational  Techniques  for  Fluid 
Dynamics,"  Vol  I,  Springer-Verlag,  New  York  1990. 


37 


15.  Sorenson,  R.L.,  "A  Computer  Program  to  Generate  Two- 
Dimensional  Grids  About  Airfoils  and  Other  Shapes  by  the 
use  of  Poison's  Equation,"  NASA  Technical  Memorandum 
81198,  1981. 

16.  Steinbrenner, J.P. ,  Chawner , J.R. , "The  Gridgen  V.8  Multiple 
Block  Grid  Generation  Software,"  MDA  Engineering,  Inc., 
Arlington,  Texas,  December  1992. 

17.  Bunning,  P.G.,  and  others,  "Overflow  User's  Manual  Version 
1.6g,"  NASA  AMES  Research  Center,  California,  June,  1992. 

18.  Shapiro,  A.H. , "The  Dynamics  and  Thermodynamics  of 
Compressible  Fluid  Flow, "  The  Ronald  Press  Company,  New 
York  1954. 


38 


w 

0) 
-H 
T! 

O 
J3 

0) 

(U 

1— 1 

a 
e 

0 

u 

5-1 
Q) 

4-) 

i-i 
o 

O 

H 

(Ti 
X 

0 

o 

00 

KD 

ro 
H 

CO 
H 

CN 

H 

ID 

H 

LD 

1 

m 

CN 

1 

CN 

r-i 

1 

r-l 

H 

H 

r-i 
r-i 
(J\ 
iH 
rH 
tH 

CTi 
CTi 

a^ 

iH 

o 

CN 

H 

rH 
H 

rH 

H 

rH 

o 

CN 
r-i 

rH 

O 
CN 

rH 

o 

CN 

rH 
O 

CN 
H 

O 
CN 
rH 

m 

QJ 

-H 

-o 
o 
£1 

0) 

JJ 

OJ 

.-1 
a 
6 
o 
u 

f-i 

o 

X 

m 

H 

X 

\> 

0 

o 

00 
H 

U) 

m 

H 

CN 

r-i 

(N 

H 

ID 

H 

LD 

H 

1 

CN 

1 

CN 

H 
1 

H 

H 
H 
H 

H 
rH 

m 

H 
H 

r-i 

r-i 

m 

r-i 

m 

r-i 

m 

H 

ro 
H 

ro 

rH 

O 
rH 
H 
H 
rH 
rH 

O 

H 

O 
VD 

O 
KO 

O 
O 

03 

-H 

0 
X! 
U 
QJ 
jJ 

o 
m 

H 

X 

ID 

0 

o 

CO 

H 

UJ 

CN 

m 

iH 
CN 

CN 

rH 

O 

un 

H 

1 

m 

CN 

CN 

r-i 
1 

H 

rH 

H 

H 

iH 

00 
LD 

r-i 

00 
ID 

00 
LD 

00 
LD 

00 

in 

00 
LD 

r-i 

iH 
H 

m 

r-i 
iH 

H 

H 

m 

rH 

m 

rH 

ro 
H 

rH 
rH 

ro 

H 

O 
VD 

rH 
H 
rH 
H 
H 

O 

rH 

o 

O 
UD 

O 
«^ 

O 

4J 
H 

Q 
H 

< 

u 

PQ 

H 

Q 

m 

H 

u 

u 

CQ 

CO 

u 

U 
PQ 

CO 

u 

CQ 

u 

CQ 

OJ 

u 

14-1 

en 

03 

o 

£1 


O   0) 

4-1     > 
-H 

>^^ 

?    0) 

05  a 
^  en 

o  '^ 
■"'  i" 

?^^ 
o 


14-1 

o 

c 
o 

-H 

JJ 
u 

0) 

•H 

(U 

4J 

03 
C 
-H 
Tl 
5h 
O 
O 
U 

II 

Pi 
H 
Q 
CQ 


4J 

C 
0) 

en 

(U 

dj 
u 

m 


CN 


ro 


OJ 


o 

-H 

4-) 

u 

QJ 
U 

-H 

T3 

QJ 
> 

-H 
4J 

03 
IT) 
QJ 
C 


O 
iJ 

QJ 
> 

•H 
4J 
-H 

cn 
o 


6 
O 

M-l 

cn 

QJ 

e 
o 
u 

o 

rH 
U-l 


O 
O 
CO 


cu 


5 

O 

rH 
14H 
4-) 

O 

o 


H  >  CO  C/] 


CUrHLDCNrOinOCNO 
>H  H   H   H  ro  ro  "t^ 

CQ 

H 


39 


03 

H 
Q 
O 
pq 
(4 

Eh 

< 

Pi 
O 

PE4 


u 

-H 
C 
O 

u 

u 

CO 

M 
> 

LT) 

in 

O 

in 

o 

H 

o 

00 

o 
o 

M 

O 
00 

H 

ID 
CO 

H 

in 

CO 
H 

O 

cn 

o 
o 

cs 

o 

CN 
CN 

rH 
1 

CO 

u 

CO 
H 

> 

O 
H 

in 

o 
in 

H 

o 

H 

in 

H 

o 
en 

H 

O 

in 

o 

00 

H 

in 

00 

H 

o 
o 

CN 

o 

CN 
CN 

H 

H 
1 
CO 
Pli 
CO 

U 
CO 
H 
> 

O 
H 

in 

O 

in 

H 

in 
in 

H 

in 

H 

LD 
00 

rH 

H 

O 

in 

o 
00 

H 

in 

00 
H 

in 
en 

H 

o 

04 

O 

H 

1 

CO 

CO 

u 

CO 

H 
> 

O 
-* 

H 

o 

o 
in 

H 

in 
in 

H 

o 

o 
00 

H 

H 

O 

O 
00 

H 

in 

00 
H 

o 
en 

H 

o 
o 

CN 

o 
m 

CN 

8 
U 

IT) 

o 

H 

O 
H 

H 

in 

H 

o 

H 

o 
m 

H 

o 
in 

H 

o 

2 

^1 
o 

u 
13 
o 

o 
u 

C! 
O 

-H 

:3 
iH 

o 

CO 

c 
o 

-H 

a 

Xi 

u 

iJ       IM 

5-1 

OJ      T3 


a* 


B 

CO 


o 


-H 
U 

m 

-H 

> 


5 

o 


13 
o 
u 
en 

■H 

> 


1 

U 

CO 

g 

CO 
H 

,, 

CO 

H 

> 

(U 

OJ 

0 

s 

40 


W 

o 
u 

1 

E-" 
< 

(D 
C 
•H 

Cm 

O 

Lf) 

o 

O 

m 

o 

in 

o 
in 

rH 

O 

rH 

O 
rH 

o 
cr> 

rH 

m 

O 
U 

O 

in 
o 

o 
m 

H 

O 

O 
LD 

rH 

in 
in 

H 

O 
rH 

O 

O 

cn 

H 

w 
o 

1 

(U 

m 

u 

0 

u 

o 
c^ 

o 

LD 

o 

o 
in 

H 

in 

U3 

rH 

in 

in 

rH 

in 

00 
H 

O 

o 

rH 
C^ 

CN 

iH 

1 

CO 

CO 

1 

c 

-H 

PL4 

O 

o 

H 
H 

o 

H 

in 

o 

H 

o 

H 

O 
00 

rH 

o 

rH 

o 

H 

CNJ 

(U 

5-1 

0 
U 

O 

LD 
O 

in 

H 

o 

H 

o 

H 

O 
CD 

H 

O 

cn 

H 

o 

H 

O 

m 

CN] 

H 

H 
1 
CO 
04 
CO 

1 

Eh 

G 
-H 
[I. 

O 

O 
H 

o 

H 

in 

rH 

o 

rH 

O 
rH 

o 

00 

rH 

o 
cn 

rH 

o 

H 
CN 

Q) 

cn 

u 

rd 
0 

u 

O 

LD 
O 

H 

in 
m 

H 

o 

KD 

H 

O 

in 

CO 

H 

in 

CTl 
rH 

O 
rH 

CN 

O 

cn 

CN 

Q 

H 

o 

s' 

O 

in 
o 

LD 
O 

in 
o 

in 
o 

in 
o 

H 

O 

H 

H 

O 
CN 

rH 

O 

m 

H 

o 
in 

H 

41 


o 


o 


i 


> 


R 


CO 


^ 


m 


en 


CN 


CO 

I 


CN 


oc 


o 


<N 


(^)IS>I 


42 


o 


o 


^ 


C/3 

> 

8 


1 

r 

—f 

o- 

O 
IL 

j  

s 

/ 

1 

>: 

— 

1 

1               i 

oo 


VO 


CN 


to 


in 


in 

CO 


m 


vn 
ri 


-      04 


C/5 


J*- 


^ 

?JI 

wn 

%-• 

'U 

CO 

& 

o 

»— 

•s 

g 

vn 

s 

o 

1-^ 

8 

"i 

o 

z 

H 

(J)B19Z 


tc 


43 


00 

w 

< 


CO 


:z; 

D 
O 


U 

o 

CO 

m 
O 

CO 


L_ 

en 


_1_ 

oo 


ri 


ri 


ri 


OJ 


oo 


r -    m 


<N 


r^ 


?x 


§ 

CO 


•c 


.1 

5 


44 


CO 

W 
U 

< 


00 

>- 

Q 

Z 
D 
O 

u 

I— t 

z 

o 

CO 

a: 
w 

CO 


I 

1  :     • 
1  :     ' 

'  : 
1  :     . 

J V-: 

IT)  ON^  — 

m  vci  CO 
ON  M3  T^ 

•          •          • 

!           !           ! 

II   II;  II 

,o  ,d  ,o 

Jl-  Ji_,  Ju- 

».         IK         ^ 

""■"vo  Cii  O 

O  — 1  <N 

IL  111  IL 

Iff 1 

'                                     *                                         « 

1                                     '                                         \ 
•                                         * 

: 

i 

i             1 

i 

i             i 

-   T^ 


m 


wn 


in 


ri 


in 


o 


(N 


IT)  -^ 


in 

o 


vn 

On 


o 


in 

oo 


OJ/i 


45 


46 


47 


48 


•5 


a 

8 


Q 

en 


ON 


49 


50 


t 


1 

8 
t 

u 

ij 

s 
8 

Q 


u 

■I 


51 


52 


53 


/       / 


/    ; 


/     / 


/       / 


o 


/  /'  /  :  ■• 


r  _'■'    /    ' 


■^  '  '■/y/-''-K'-A/i 


i^    ■•'     /■    /     /    /-   /'   //    //,-^A..     '•■ 


//         /  r       '        //  f       *■'/'•      f     -^     ■'     /  , 


1 1 


\   o 


'      A      -    '^  ;,         A^^   A.^'^  A  xr>^y^/^AA 

•        ,  -"       A.       X  ^-    /^    /^    /^    /■   /■   '     / 

-.  V        •  ■  A  -^^^-      ^     X    /<     ,-      A\//    /-    / 

^       •■  .--       .-.      ,       .'•-     ,-.     A-     /     /     /    /    /    / 


/-•      /^ 


U 


54 


55 


CO 


a 


o 

g 
i-i 


00 

E 


56 


2<cS 


W  Q 

^2 


O  VO  (N 

Q  o  o  o 
o  -^  -^  o 


U 


o 

1— ( 

II 

JZ 

o 

s 

o 
c: 


•a 
O 


2 


tt^^,^-J^i 


U] 

> 

3 

OS 

D 

™ 

";-:; 

w 

c^ 

~ 

C;;, 

— 

C; 

c 

cs 

X' 

z, 

C  O 

f_,0 

;3r:^.  5s: 

sr 

»K 

=r. 

■w 

■s: 

.T< 

sr. 

■at 

•x 

.afs:; 

^  ^ 

•-/'*'. 

Ci/' 

CT 

^.z- 

*•* 

•r. 

^ 

■<'. 

~ 

<r-. 

'^^' 

\r. 

3'^- 

OH 

cr 

— '  — 

c . 

'■• 

r*- 

r^. 

**?■ 

-r 

k/-- 

V; 

w 

•^ 

r-  r- 

UO 

— 

OCT 

■z. 

-- 

- 

^ 

c; 

-^ 

— 

^ 

c- 

- 

__ 

>■■:!/'/ 


/y"\/ 


sSBSSBSS 

O  i/^C  IV,  O  i/~.  O  <r, 
<i-,  lA:  vo  o  r- r- oc  00 


57 


58 


CSJ 


o 


On 

E 


59 


5S 


OSpO 


o  ^  <~< 

Qooo 

.  00  H  9  2 


S 


II 

s: 
o 

5 


5  -2 


o 

o 

c 


■n 
O 


2 

a 


._.j 


•9 

cq 


O 

B 


10  0 10  Q  ir,  o  ir)  O 
■^  m  10  vo  vo  r~- r- X 


O 


60 


5^ 


<u  £  PS 


O  ^  r-i 

wi  : 

Q  o  o  o 

^  00  H  H  2 
-H  -I  vo  -:f  .-r 


(S 


5  ^ 
z  a 

U    0 


■o 

■c 
O 


_L. 


r^x 


C3 


C/3 

o 

a 


oooo 


in  Q '/^  C  lo  O  (/",  c^ 
.'o^-or-r-ocoocF. 


(D 


E 


61 


to 


5=) 

8 


C/3 

e/3 

IS 

O 

a 


62 


o 


C/3 

-a 

o 

g 


CM 

oo 


63 


I 


'^    ^ 


2  > 


—  < 


n  ^o  CM 

o  S  ^  ^  ^^ 
"~.  Sc  ^.  ^  -jc 


1/3 

-} 

U 

> 

u 

-] 

/v* 

3 
OSS 

;.~ 

';^ 

^^. 

~  ~ 

•JiK. 

- 

^' 

Cl-r 

^jk: 

-~.5  ;■ 

Ho^ 

"«. 

;Z^ 

."Z.' 

.X.'  va\ 

TX^. 

.■.-' 

.'— ; 

'«* 

^jsr,  rjsr. 

,'— . 

^^  ^% 

ZO^r 

:~ 

V", 

-'*.'w 

v^ 

■/ 

■■«* 

=  i'^  vJ" 

Z'AC;   ": 

— 1 

— 

•■■■■ 

^•■■^r*'" 

•"" 

""' 

/  _ 

v^:-C 

V, 

■■~~v~- 

'- 

— 

■-" 

" 

■"■ 

"■ 

.^. 

■w 

_ 

_ 

... 

0  3 

cr 

u-.O'-r.c: 

;c-. 

OV, 

'r 

in 

.^ 

Z;_ 

— 

iCOC 

-6 


•G 


o 


1 

■a 

5 


CO 

o 


64 


_ 

< 

< 

< 

(U 

2i 

HO 

n 

vO 

n 

* 
* 

* 

^ 

^ 

"n 

f^ 

:'  : 

X 

X 

^ 

'    . 

oo 

■*— 1 

^-t 

1/^ 

— 

— 

vC 

r-i 

\ 


\ 


11 


<     o 


> 


\ 


OS: 


O  V-,  c  •  r ,  w 'vn  o 


•8 

CO 


O 


I 

Si 


65 


ijcif  Jlt 


:r-t  (N  (N  VO 


I    I    I    loio 


lA-.Ui  CO  CO 


H 
0) 

•o 

0 

B 

>i 
TJ 
0 
i) 
M 
4) 
+J 
«W 
(0 

0 


u 


u 
o 

I 

o 


o  >  o 

00   B   (Q 
•H  H   H 


;  /( 


'\  III 

; 

:      III 

:              / 

i  //// 

:     'it* 

i            / 

Mil! 

i          1 

Mil:  I 

■.ill   '    1 

'.         1 

]■    f 

:     / 

m 

•    1 

i  /             •• 

m 

A     lit 

'I 

.'  \       //• 

/          //.' 

I: 

^ vm 

*  \ 

j  /  //  / 

• 

f  /•'  / 

! 

;/  //  / 

: 

Mr 

i 

//'V 

: 

ft, 


m 


n 


O 


<N 


//  ;.• 


M 


€1 


<^  d' 


o 

I 


fN 


O 
I 


o 
I 


00 

o 

o 
I 


o 


o 
I 


o 


o 
I 


o 


o 
I 


sssa<3-po 


66 


1 

^p^ir 

•        B  ^ 

/ 

jj*  Oh       ^  R 

'^I^^^^Hk^ 

/ 

S  i-J  a>  H  M 

^^^^^^bL 

S<p!;Po     jv^i^^ 

ilJj^^^^^^^E' 

/•■ 

O^Cs> 

H^^K 

///// 

^^1 

^^^^HHIHI*' 

Q  o  o  o                               .... 

1      1  / 

OO   X   X-ri                                                             / 

III 

'^0°'^<S                                                         ,  :  ■-    .            / 

////'. 

If, 

^^^  - 

^""^-^...^"  y 

m 

1       1       1 

1 

N 1 — ■■ ' 1^, -J^ 

rr^^sB 

y             '^^  /'I'i 

{      (  'v 

■■mjhbmH 

■■■i^ihB 

^ 

■■^■i^H 

lo 

V 

vrnmi.,:^^ 

00 

^*\ 

\ 

■■^■^■■1 

'-gWBUR^-*'^-^-- 

^ 

_fl 

\ 

, 

■SBBBI 

rf'i^^^^^ 

"S 

'v 

, 

■■^■laaH 

cd 

•v. 

T 

SOBRtt^;./.  V. 

^^ 

\.                                       --..^ 

J 

.IEMSScf!:a^\ 

mtj 

s 

- 

iSSSSI 

IP^^'IH 

1 

ssssi 

■_  ^^ 

UiZM 

^1 

■ .  -  jriw^^^BSnK  - 

1^ 

4^ 

whihi^hh 

^^^H^N* 

a 

; 

ih^b^h^hH 

'^HRs^ 

D 

1 

■BSSSi 

^HnL, 

Z 

t+in 

■^B^^^mI 

hri 

A 

^^Ki'sT'H 

^^ 

Kj 

T3 

: 

"■  — 

^^■^^(^fl 

hkI 

O 

IS 

WBKKt- 

El 

^ 

^ 

.^raKS» 

^ 

'"  ^'MSHBH' 

Q> 

^ 

:•;<.^M<3^B■ 

(L> 

^ 

^^mtf.KF'-          -i^m 

_ 

^ 

V^hS^hCuI^h^  • 

*&. 

^^^''^ 

: 

H^^'Vs^^Jl 

a 

^' 

o 

X 

■ 

iiiii 

' '^*v*H*H^^^- ' 

o 

11 
a> 

•  {   ( 

^^^^^. 

CO 

\                  \ 

•^-    ■ 

3 

o 

\    \ 

m 

^^^^^^^^^K     '" 

.''^H 

^v     ^ 

^^^H 

:^'2^ 

''^^^1 

^^u 

_,_———""■   '      '"~ 

^^^^Bk^^WBH^ 

-^^^^^^^B 

^■K.%*mh 

^^^&>\*«SS 

^^■^'^«»bS 

v 

I 

t 
< 

> 

ooooo 

1 

1 

-  -  rrcooooooo 

oooooooooooooooo 

OOOOOOOOQOOOOOOO 

'^1 

^K 

'  — —  oooocoooo 

(^MlfiiJl»ayiaTjiwirail^UJil^Mr^T»l 

c 

i~  icro'ooocJooooo'oocSocJOOT-iT-J^T-^T-H^T-l 

''IfflH 

l^p 

Vm 

00 

o 

«5 


C/3 
Oh 
C/3 


T3 

O 

u 
B 


u 

*-> 

« 

3 

a 

CO 

u 

3 


69 


70 


o 


o 


Oh 
CO 

I 

T3 
O 
X) 

<u 

E 


ed 

3 

o 

ri 

m 

u 


CN 


o 

S 


u 

4J 
C 
U 

x» 
u 

o 
o 

(0 

-H 


3 
XJ 
-H 
(0 
Q) 
PC 


o 

I 

9> 


4-)  4J  4J 

3  3  3 

O  O  O 

rH  CM  n 

tT>  tT>  tj> 

W  (0  M 

M  M  M 


O 
O 


O 
O 
ID 


O 
O 


CO    (U 
4J 


O 
O 
CNJ 


O 
O 


O 
I 

0) 


4= 


CO 

CO 

t 


o 

c^ 

o 
•a 

pes 
u 


73 


o 


o 


3 

a 


"" 

"! 

'          1 

1          i                   •              1    >         <                   •              1    '          •                   '              J    ' 

■  / 

1 

! 

/' 

1 

/'' 

-  ♦ 
( 

J  4J 
3   3 
3   0 

*          » 

♦J 

9 
0 

/' 

<N  •*  VO 

y 

iH  fo  in 

/ 

' 

•M  M  <M 

/ 

-      C 

Q    CO 

a 

/ 

- 

1 

i 

4   U 

M 

*• 

- 

* 

/ 

o 

/ 

«-• 

• 

/ 

ri 

ja 

/ 
/ 
/ 

u 

/ 

13 

I 

/ 

/ 
/ 

/ 

g 

/ 

/ 

O 

/ 

M 

/ 

Xi 

/ 

U 

/ 

/ 
/ 

«H 

— 

• 
/ 

— 

0 

/ 

/ 
/ 

tr 

/ 
/ 

0 

/ 

p 

/ 

CO 

_ 

/ 

- 

•H 

/ 

n 

/ 

H 

/ 

(0 

/ 

3 

/ 

-0 

_ 

_ 

•H 

la 

0 

(x: 

1 

/ 
\ 

'    1          " 

■  7^^    r~  III            1            — u_ 

-1 1 — 

o 
I 


o 
o 
in 


o 
m 


o 

o 


o 
in 


o 
o 

00 

c 

0 

O    10 

in  ki 

(N   0) 


o 
o 


o 
in 


o 
o 


o 
m 


1»= 


00 

o 
I 


o 


-a 

•a 

? 
t 

u 
u 

s 
8 

o 
o 


C3 

u 


75 


< 

i      ; 

.  1  m 

1                   1 

•^    .-)    CNJ 

t 

rH         .         • 

t 

_..  .  .•   t-H  >H 

; 



U       1      1 

i 

^^.^ 

1    CO  to 

! 

T3 

r-i     On     CU 

• 

-H 

ni  m  in 

I 

o 

H 

c 

f!?+ 

0) 

o 

/ '.  ' 

c 

o 

'.1 

•H 

1  ••  ' 

<M 

1  •  *j 

— 

/  '*  V 

~ 

>, 

//I 

X) 

/<  ,'1                   ■ 

0 

it.'  • 

XI 

T^i 

0) 

1 '  • 

»-) 

/  1  '■ 

(U 

1  1  • 

rH 

Aj-l 

a 

*" 

g 

^y-  '   • 

o 

^^  .'  ''    '• 

o 

y^       :   1        \ 

u 

< 

ri-/        1 

o 

^^ 

<M 

^^^^^^%4.^ 

^t 

^^"^^***--^  *  •"*»■ 

JZ 

u 

• 

■'*'^-^-*w 

(0 

s 

: 

■--=^;a 

CO 

> 

i 

"\N 

M 

: 

\    \ 

(0 

i 

\    \      - 

>-l 

Q. 

; 

\  > 

-d 

i 

\ 

o 

j 

t 

\ 

\ 
1 

i 

\ 

1 > « 

in 
o 


o 

(0 


oo 
O 


<3 


T 
O 


o 


CM 

o 


o 
o 


ss3Jd- p3 


76 


APPENDIX 
COMPUTER  PROGRAMS 
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PROGRAM  KSI.F 
*********************************************************** 

*  This  program  is  developed  to  calculate  KSI (X) ,  K(X)  and  X  * 

*  using  mumerical  integration  based  on  trapezoidal  rule  to  * 

*  solve  equation  (32)  and  plotting  the  output  as  shown  in  * 

*  Figs .  1  .  * 
************************************************************* 

REALMO),  X(0:401,3),  A,  Al ,  B(401,3) 
#     ,  KSI,  P,  FUNG,  H,  DD,  Q,  L,  K(401,3),  N,  Nl 
INTEGER  YY 

OPEN (UNIT=9 , FILE= ' KSI ' , STATUS= ' UNKNOWN' ) 

PRINT  *,  'ENTER  LOWER  &  UPPER  BOUND,  #  OF  INTERVALS,  #  OF 
#DATA  SET' 
READ  *,  A,  Al,  N,  Nl 
100   DO  10  1=1,3 

PRINT* ,  ' ENTER  MACH  NO . ' 
READ*,  M(I) 
C  M(I)  =  0.9+1*0.1 

IF  (Md)  .LT.1.0)  THEN 
P  =  -1. 
ELSE 
P  =  1. 
ENDIF 
DO  20  J=1,N1 

B(J,I)  =  J*A1/N1 
H  =  (B(J,I) -A)/N 
AREA  =  0. 

K(J,I)=(-0.0102/(M(I) **2) ) * ( (ABS( (B(J,I) **2) - ( (ABS(1- 
#   M(I) **2) ) **1.7574) ) ) ** {2./3.) +(ABS(1-M{I) **2) ) **1.1716) 
C 
,  DO  3  0  L=  1,N 

KSI  =  A  +  (L-0.5) *H 

DD  =  (KSI**2+P* (ABS (l-M(I) **2) ) **1.7574) 
IF  (DD.GT.0.0)  THEN 
FUNG  =  1./ (DD** (1./3.  )  ) 
Q  =  1. 
ELSE 

FUNG  =  1./ (  (-l.*DD) **  (1./3. )  ) 
Q  =  -1. 
ENDIF 

AREA  =  AREA  +  H*FUNG*Q 
X(J,I)  =  AREA 
3  0       GONTINUE 
2  0       GONTINUE 
10       GONTINUE 


C 


DO  4  0  J=  1,N1 

PRINT  50,   (B(J,I) ,X(J,I) ,K(J,I) ,1=1,3) 
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50       F0RMAT(1X,3  (F6.3,2F10.5,1X)  ) 
40    CONTINUE 
c 
2  00   PRINT*, 'TRY  AGAIN?,  ENTER  1  FOR  YES , 2  FOR  DATA  FILE, 
#  OTHERS  FOR  NO' 
READ*,  YY 
IF  (YY.EQ.l)  THEN 

GO  TO  100 
ELSE 

IF  (YY.EQ.2)  THEN 
DO  110  J=  1,N1 

WRITE  (9,50)  (B(J,I) ,X(J,I) ,K{J,I) ,1=1,3) 
110  CONTINUE 

GO  TO  200 
ENDIF 
ENDIF 
END 


PROGRAM  ZETA.F 
************************************************************ 

*  This  program  is  written  to  calculate  ZETA(r)  and  r  using* 

*  equation  (29)  and  plotting  output  in  Fig. 3  .  * 
************************************************************ 

REAL  M ( 3 ) ,  ZETA (0:14,3),  R ( 14 ) 
0PEN(UNIT=9,FILE='ZZ' , STATUS =' UNKNOWN' ) 
C 

DO  10  1=  1,3 
C        M(I)  =  0.7+1*0.1 

PRINT* ,  ' ENTER  MACH  NO . ' 
READ*,  M(I) 
DO  2  0  J=  1,14 
ZETA(0,I)  =  1000.0 
PRINT*,  'ENTER  R  VALUES' 
READ*,  R(J) 
C         R(J)  =  0.1*J 

ZETA(J,I)  =  (4/(R(J)**2) )+(ABS(l- (M(I) **2) ) 

#  *(R(J) **2.8284) ) 

#  +( (l-M(I) **2) * (R(J)**7.657)/50.63 

IF  (ZETA(J,I) .GE.ZETA(J-1, I)  ZETA ( J, I) =ZETA ( J- 1 , I) 
2  0      CONTINUE 
10    CONTINUE 

PRINT  30,  M 
30    FORMAT (11X,5F10. 4) 

DO  40  J=  1,14 

WRITE(9,50)   (R(J),   (ZETA (J, I) , 1=1 , 4 ) 
50    FORMAT(4X,3 (2F10.4,2X) 
40    CONTINUE 

END 
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PROGRAM  D2D3 . F 
*********************************************************** 

*  This  program  will  read  data  files  output  Gridgen2d    * 

*  from  2D  to  3D  for  further  processing  with  rotategr.f  * 
*********************************************************** 


REAL  X(77,60) ,Y(77,60) ,Z(77,60) 
READ (12.*)  IDIM,JDIM 


# 


READ (12, *) 

Z(I, J) =0.0 
WRITE (14,* 
WRITE (14, * 


(X ( I , J)  , 1  =  1 , IDIM)  , J=l , JDIM)  , 
(Yd,  J)  ,1  =  1,  IDIM)  ,J=1,JDIM) 

IDIM,  JDIM,  1 

( (X(I, J) ,1=1, IDIM) , J=l, JDIM) , 
( (0.0, 1=1, IDIM) , J=l, JDIM) , 
(  (Yd,  J)  ,  1  =  1,  IDIM)  ,  J=l,  JDIM) 


STOP 
END 


PROGRAM  HALF.F 
************************************************************* 

*  This  program  will  write  out  half  of  a  grid  file  output   * 

*  from  Grape  for  further  processing  with  rotategr.f       * 
************************************************************* 


DIMENSION  X(200,200) , Y(200,200) 
READ (30,*)  IDIM, JDIM 
READ  (30,  *)   (  (Xd,  J)  ,1  =  1,  IDIM)  ,  J=l, 
#  ( (Y(I,J) ,I=1,IDIM) ,J=1, 


JDIM) 
JDIM) 


IDMD2=IDIM/2  +  1 


10 


DO  10  J=1,JDIM 

Y(IDIM2,J)  =0.0 
CONTINUE 


WRITE (31, 

* 

WRITE (31, 

* 

# 

# 

STOP 

END 

IDMD2, JDIM,1 

( (X(I, J) ,I=IDMD2) , J=l, JDIM) 
( (0.0,I=IDMD2) , J=l, JDIM) , 
(  (Yd,  J)  ,I  =  IDMD2)  ,  J=l,  JDIM) 
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PROGRAM  ROTATEGR.F 
************************************************************* 

*  This  program  will  create  3D  grid  (body  of  revolution)    * 

*  by  rotating  3D  plane  grid  for  further  processing  with   * 

*  Flow  Solver  * 
************************************************************* 

DIMENSION  X(200, 200) ,  Y(200,200),  Z(200,200) 
DIMENSION  XX(200, 100,200) ,  YY (200 , 100 , 200 ) , 
#ZZ(200, 100,200) 
CHARACTER* 3  0  FNI 
CHARACTER* 3  0  FNO 


PRINT* 


READ(* 
PRINT* 
READ(* 
PRINT* 
READ(* 


' INPUT  GRID ' 


21)  FNI 

'SG  =  0,  MG  =  1',  ',UNF=1,  F0RM=2' 
* )  IGRI 

'IFORI  =  1  UNF,  IFORI  =  2  FORM' 
*)  IFORI 
21    FORMAT  (A) 

IF  ( IFORI. EQ.l)  THEN 

REWIND  1 

OPEN  ( 1 , FILE= FNI , FORM= ' UNFORMATTED ' ) 

IF  (IGRI. EQ.l)  READ  (1)  MGR 

READ  (1)  I1,J1,K1 

READ  (1)  ( (X(I, J) ,1=1,11) ,J=1,J1) , 

#  ( (Y(I,J) ,1=1,11) ,J=1,J1) , 

#  (  (Z(I,J)  ,1  =  1,11)  ,J=1,J1) 
CLOSE  (1) 

ENDIF 

IF  (IF0RI.EQ.2)  THEN 

REWIND  2 

OPEN  (2, FILE= FNI, FORM=' FORMATTED' ) 

IF  (IGRI. EQ.l)  READ (2,*)  MGR 

READ (2,*)  I1,J1,K1 

READ (2,*)  ( (X(I,J) ,1=1,11) , J=1,J1) , 

#  ((Yd, J)  ,1  =  1,11)  ,J=1,J1)  , 

#  ((Z(I,J) ,1=1,11) ,J=1,J1) 
CLOSE  (2) 

ENDIF 

PI  =  4.*ATAN(1.) 

PRINT*,  'NO  OF  PLANE  IN  J  DIRECTION  ?' 

READ ( 5 , * )   JM 

DTH  =  (180./ (JM-1) )* (PI/180) 

DO  11  1=1,11 

DO  11  J=1,J1 

K=J 

XX(I,2,K)  =  X(I,K) 

YY(I,2,K)  =  0.0 

ZZ(I,2,K)  =  Z(I,K) 


11   CONTINUE 
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20 


30 


IM=I1 

KM=J1 

DO  20  J=3,JM+1 

DO  20  1=1, IM 

DO  20  K=1,KM 

XX(I, J,K) =X(I,K) 

TH  =  (J-1) *DTH 

YY(I,J,K)  =  SIN(TH) *Z(I,K) 


ZZ(I, J,K) 

=  COS (TH) *Z(I 

CONTINUE 

DO  30  1=1 

IM 

DO  30  K=l 

KM 

J=l 

XX(I, J,K) 

=  XX(I,J+2,K) 

YY(I,J,K) 

=  -YY(I,  J+2,K) 

ZZ(I, J,K) 

=  ZZ(I,J+2,K) 

J=JM+2 

XX(I, J,K) 

=  XX(I,J-2,K) 

YY(I,J,K) 

=-YY(I, J-2,K) 

ZZ(I, J,K) 

=  ZZ(I,J-2,K) 

CONTINUE 

JM=JM+2 

PRINT*,  'Output  filename  =' 
READ (5,21)  FNO 

REWIND  3 

OPEN(  3,  FILE=FNO,  FORM= ' UNFORMATTED ' ) 

WRITE (3)  IM,JM,KM 

WRITE(3)  (((X(I,J,K),  I=1,IM) , J=l, JM) ,K=1,KM) , 

#  (((Y(I,J,K),  I=1,IM) , J=l, JM) ,K=1,KM) , 

#  (((Z(I,J,K),  I=1,IM) , J=1,JM) ,K=1,KM) 
CLOSE  (3) 

STOP 
END 
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PROGRAM  READX.F 
************************************************************ 

*  This  program  will  read  3D  grid  file  and  convert  from   * 

*  formatted  into  unformatted  file  (grid. for  to  grid. in)  * 

*  for  further  processing  with  Overflow.  * 
**************************************************************** 

DIMENSION  X( 77, 13,60) ,  Y(77,13,60),  Z(77,13,60) 
OPEN (UNIT=12,FILE=' GRID. FOR' , STATUS =' UNKNOWN' ) 
0PEN(UNIT=14,FILE='GRID.IN' ,STATUS='NEW' , 
#FORM= ' UNFORMATTED ' ) 
READ (12,*)  IDIM, JDIM,KDIM 
READ(12,20)  (((X(I,J,K),  1=1 , IDIM) , J=l, JDIM) , K=l , KDIM) , 

#  (((Y(I,J,K),  1=1, IDIM) ,J=1, JDIM) ,K=1, KDIM) , 

#  (((Z(I,J,K),  1=1, IDIM) ,J=1, JDIM) ,K=1, KDIM) 
c 

WRITE (14)  IDIM, JDIM, KDIM 

WRITE(14)  (((X(I,J,K),  1  =  1, IDIM)  ,J=1,  JDIM)  ,K=1, KDIM)  , 

#  (((Y(I,J,K),  1=1, IDIM) ,J=1, JDIM) ,K=1, KDIM) , 

#  (((Z(I,J,K),  1=1, IDIM) ,J=1, JDIM) ,K=1, KDIM) 
20    FORMAT (6F15. 6) 

STOP 
END 


PROGRAM  READQ.F 
************************************************************ 

*  This  program  will  read  q.save  file  output  from       * 

*  Overflow  and  convert  from  unformatted  into  formatted  * 

*  file  for  further  processing  with  Plot3d.  * 
************************************************************ 

DIMENSION  0(77,13,60,5) 

OPEN (UNIT=1,FILE=' Q.SAVE' , STATUS= ' OLD ' , 
#FORM= ' UNFORMATTED ' ) 
OPEN(2  0,FILE='Q.FORM' , STATUS= ' NEW' , FORM= ' FORMATTED '  ) 
READ(l)  NI,NJ,NK 
READ(l)  FSMACH,ALPHA,RE,TIME 

READ(l)  ( ( ( (Q(I,J,K,NX)  , 1=1, NI) , J=1,NJ) , K=1,NK) ,NX=1, 5) 
C 

WRITE(20,*)  NI,NJ,NK 

WRITE (20,*)  FSMACH,ALPHA,RE,TIME 

WRITE(20,*) ( ( ( (Q(I, J,K,NX) 

#,I=1,NI) ,J=1,NJ) ,K=1,NK) ,NX=1,5) 
PRINT*,  NI,NJ,NK 
STOP 
END 
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PROGRAM  OVERFLOW. IN 
************************************************************ 

*  Input  file  for  Overflow  to  run  either  EULER  or  NAVIER-   * 

*  STOKES  Flow- solver.  * 
************************************************************ 

$GLOBAL 

CHIMRA=.F.,  NSTEP=100,  RESTART=.F.,  NSAVE=100,  NQT=0, 
$END 
$FLOINP 

ALPHA=180.0,  FSMACH=0.95,  REY=  6.00E6,  TINF=  520.0, 
$END 


77x13x60 


$GRDNAM 

NAME=  'AXI- SYMMETRIC  BODY 

WITH  AFT 

-CONE,  77 

GRID' , 

$END 

$METPRM 

IRHS  =  0 

,  ILHS  = 

=  2, 

IDISS 

=  2, 

$END 

$TIMACU 

DT  =  0.1 

,  ITIME 

=  1, 

TFOSO  =  1.00, 

CFLMIN  = 

$END 

$SMOACU 

ISPECJ  = 

2,  DIS2J  = 

2.00, 

DIS4J  = 

0.02, 

ISPECK  = 

2,  DIS2K  = 

2.00, 

DIS4K  = 

0.02, 

ISPECL  = 

2,  DIS2L  = 

2.00, 

DIS4L  = 

0.02, 

SMOO  =  1 

.00,  EPSE  = 

0.35, 

$END 

$VISINP 

VISCJ  = 

.F.  ,  VISCK 

=  .F. 

,  VISCL 

=  .F.  , 

NTQRB  = 

0, 

ITTYP  = 

1, 

ITUIR  = 

3, 

JTLS   = 

1, 

JTLE   = 

77, 

KTLS   = 

1, 

KTLE   = 

13, 

LTLS   = 

1, 

LTLE   = 

60, 

TLPAR1= 

.3, 

$END 

$BCINP 

NBC 

6, 

IBTYP  = 

15, 

12, 

12, 

1,   32, 

15, 

IBDIR  = 

1, 

2, 

-2, 

3,    3, 

-1, 

JBCS   = 

1, 

1, 

1, 

1,    1, 

77, 

JBCE   = 

1, 

77, 

77, 

77,   77, 

77, 

KBCS   = 

1, 

1, 

13, 

1,    1, 

1, 

KBCE   = 

13, 

1, 

13, 

13,   13, 

13, 

LBCS   = 

1, 

1, 

1, 

1,   60, 

1, 

LBCE   = 

60, 

60, 

60, 

1,   60, 

60, 

$END 

5.00, 
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